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ABSTRACT 

pi ^ Wc present 1.2 mm MAMBO-2 observations of a field which is over-dense in Lyman 

I Break Galaxies (LBGs) at 2; ~ 5. The field includes seven spectroscopically-confirmed 

O LBGs contained within a narrow {z = 4.95 ± 0.08) redshift range and an eighth at 

z = 5.2. We do not detect any individual source to a limit of 1.6 mJy/beam (2 x rms). 
^ When stacking the flux from the positions of all eight galaxies, wc obtain a limit to 

1 I the average 1.2 mm flux of these sources of 0.6mJy/beam. This limit is consistent 

with FIR imaging in other fields which are over-dense in UV-bright galaxies at z ~ 5. 

■'"^ Independently and combined, these limits constrain the FIR luminosity (8-1000/im) to 

^ a typical 2; ~ 5 LEG of LpiR < 3 x 10"^^ Lq, implying a dust mass of Mdust ^ 10^ Mq 

(both assuming a grey body at 30K) . This LpiR limit is an order of magnitude fainter 

than the LpiR of lower redshift sub-mm sources (z ~ 1 — 3). We see no emission from 

any other sources within the field at the above level. While this is not unexpected 

given millimetre source counts, the clustered LBGs trace significantly over-dense large 

scale structure in the field at 2; = 4.95. The lack of any such detection in either this 

or the previous work, implies that massive, obscured star-forming galaxies may not 

always trace the same structures as over-densities of LBGs, at least on the length scale 

. . probed here. We briefiy discuss the implications of these results for future observations 

> with ALMA. 

• ^ 
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1 INTRODUCTION 

Lyman Break Galaxies (LBGs) form a substantial fraction 
of the observed z ~ 5 galaxy population (e.g Vanzella et al. 
2009; Douglas et al. 2009, 2010), with a comoving number 
density of 4>* ~ 10"^ Mpc"^ (at 2; ~ 4 - 6, Bouwens 
et al. 2007). They are identified via their bright UV contin- 
uum emission arising from hot, young stars in unobscurcd 
starburst regions. Considerable work has been carried out 
to investigate LBG properties which can be discerned from 
their rest- frame UV spectra and UV/optical spectral energy 
distributions (e.g Verma et al. 2007; Stark et al. 2009). How- 
ever, until now little work has been carried out to explore 
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their cooler dust and interstellar gas components necessary 
for a clearer and fuller picture of these important probes 
of early galaxy evolution. In order to fully understand star- 
formation activity at high redshift and the subsequent evo- 
lution of early star-forming galaxies, we need to observe 
their complete baryonic budget. By comparing the stellar, 
molecular gas and dust fractions of these systems we can in- 
fer their star-formation history and potential fate, therefore 
understanding the importance of this population to galaxy 
evolution in general. 

Until recently, studies of the dust content of the highest 
redshift galaxies {z ~ 5) have been limited to massive and 
rare systems in comparison to more typical and numerous 
star-forming galaxies at 2: ~ 5. These studies have primarily 
targeted quasar host galaxies (e.g Wang et al. 2008; Carilli 
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et al. 2010) and the highest redshift luminous submillimetre 
galaxies (SMGs, e.g Coppin et al. 2010; Riechers et al. 2010), 
a small subset of the SMG population which is dominated 
by galajcies at 1 < 2 < 3 (Blain et al. 2002; Chapman et al. 
2005; Small ct al. 2002). Dust emission has been observed 
in more typical gala;xies at significant redshifts (2; ~ 1 — 3) 
but usually only in lensed objects (e.g Baker et al. 2001; 
Negrello et al. 2010; Conley et al. 2011), or from rare LBGs 
with exceptionally high star-formation rates (Chapman & 
Casey 2009). Given their comparatively high source density, 
2 ~ 5 LBGs are likely to be better tracers of more typical 
high redshift star-formation than such rare and highly lumi- 
nous sources, even given the star-formation rates of 10-100 
M0yr~^ determined from the UV emission of ~L* LBGs 
{e.g. Verma et al. 2007). 

Given that these star-formation rates are typically an 
order of magnitude less than those of high redshift SMGs, 
and that the surface densities of SMGs are lower than those 
of « ~ 5 LBGs, unless there is substantial obscured star- 
formation associated with typical LBGs, they will not be de- 
tected at the levels usually reached by submillimeter surveys 
(several mjy at 850/xm- 1.2mm). To detect or constrain dust 
emission from z ~ 5 LBGs, the average source needs to be 
probed to sub-mjy levels. In addition, a large-enough num- 
ber of these sources need to be sampled to search for any rare 
LBG that has significant sub-mm flux. Until recently there 
were no published predictions of the expected sub-mm flux 
of distant LBGs based on any self-consistent modelling of the 
star-formation within the population. Recently Gonzalez et 
al. (2011) derived predicted submm fluxes for high redshift 
LBGs based upon the GALFORM semi-analytic model (see 
Cole et al. 2000). These predicted fluxes are approaching 
the observational limits of Stanway et al. (2010) (see later), 
indicating that such observations are starting to challenge 
theoretical predictions for early galaxy evolution. 

While studies of typical unlensed star-forming galaxies 
at lower redshifts {z < 3 LBGs) have yet to yield an individ- 
ual detection in the FIR, we cannot necessarily extrapolate 
these results to infer the properties of the z ^ 5 population. 
Many of the properties of the z ~ 3 sources are different 
from those of typical LBGs at z ^ 5 (Verma ct al. 2007). 
The higher redshift LBGs are also likely to be the progen- 
itors of different current-day galaxy populations than the 
2 < 3 LBGs. While the bulk of stars in present-day galax- 
ies are formed at 2 < 3 (e.g. Sobral et al. 2012), ellipticals 
found in the cores of galaxy clusters appear to rapidly form 
a significant fraction of their stars at 2 > 3 {e.g. Thomas 
et al. 2010). Consequently, the mm/submm characteristics 
of the 2 ~ 5 LBGs could differ from those of 2 < 3 LBGs, 
requiring that they are directly measured rather than ex- 
trapolated from the existing work at lower redshifts {e.g. 
Webb et al. 2003). 

In order to further explore the 2 ~ 5 LBG population 
in the mm-submm regime, wo present the results of 1.2mm 
MAMBO-2 imaging of a field selected for a high density 
of spectroscopically-confirmcd LBGs which trace an appar- 
ently highly over-dense structure in the young Universe. The 
field is drawn from the ESO Remote Galaxy Survey (ERGS, 
Douglas ct al. 2007, 2009, 2010), which identified a sample 
of 70 spectroscopically confirmed 2 ~ 5 LBGs over 10 widely 
spaced ~ 40 arcmin^ fields. 

Two of these fields display highly significant 3- 



dimensional over-densities of spectroscopically-confirmed 
UV-bright sources over narrow redshift ranges {Az ~ 0.1). 
The two structures occupy no more than 4 per cent of the 
effective survey volume, while containing ~ 30 per cent of 
the spectroscopically-confirmed LBGs. Hence, the LBGs in 
these fields are likely to trace significant large scale struc- 
tures in an early stage of their evolution, and may themselves 
contain other sub-mm luminous galaxies {e.g. Capak et al. 
2011). 

The redshift and spatial clustering of sources in these 
fields fortuitously make them ideal testbeds for simultane- 
ous observations of many LBGs and additionally probe sig- 
nificant early mass over-densities. The first of these fields 
was studied by Stanway et al. (2010) at 870pm using the 
Large Apex BOlometer CAmera (LABOCA) resulting in a 
limit of S'sTO/^m < 0.85mJy for the average emission from 
nine 2 ~ 5 LBGs. The MAMBO-2 observations presented 
here are of the second highly-over-denso field in ERGS which 
cover a further eight 2 = 5 LBGs and the over-dense struc- 
ture that seven of these inhabit. By doubling the sample 
size of spectroscopically-confirmed LBGs targeted for their 
mm/sub-mm dust emission, using a field spatially unrelated 
from the first, we aim to both confirm the limit on the omis- 
sion from a typical 2 ~ 5 LBG and to search for or limit the 
number of potential mm/sub-mm luminous outliers in the 
population. 

These observations also allow the exploration of the 
nature of 2 ~ 5 LBGs beyond what is possible with the 
standard optical/near-IR analysis (e.g. Douglas et al. 2010). 
The galaxies identified through the Lyman Break technique 
are selected for their rest-frame (presumably unobscured) 
UV emission which itself implies a minimum star-formation 
rate (and therefore, via a set of assumptions, a minimum 
FIR luminosity). However, it is entirely possible that some 
of these either contain obscured stellar populations or have 
neighbours/companions that have similarly obscured stellar 
populations (which may be more significant that the UV- 
emitting population). These observations directly target the 
source of any obscuration independent of the UV emission. 

Throughout this paper all optical magnitudes are 
quoted in the AB system (Oke & Gunn 1983), and the cos- 
mology used is 100/i///o = 70kms"^ Mpc"\ f2A = 0.7 and 

nM=o.3. 



2 OBSERVATIONS 

Observations were undertaken using the MAMBO-2 array 
(Kreysa et al. 1998) on the IRAM 30m telescope in the 
Spanish Sierra Nevada. MAMBO-2 is 117 channel bolome- 
ter with individual channels arranged in a circular pattern. 
It has a bandwidth centre of ~250GHz (1.2 mm) and half 
power spectral range of 210 to 290 GHz. The effective beam 
size is 10.7" (~ 70kpc at 2 ~ 4.95) and it has an under- 
sampled field of view of 4'. Field J1054. 7-1245 was observed 
on two separate occasions using On-the-Fly mapping for a 
total integration time of 2.2 hrs on 2009 March-15/16th 
and 5.5 hours during 2011 January-March. Each map was 
comprised of 18 azimuthal subscans of 66 s each, scanning at 
5"s~'^, with 8" steps in elevation between subscans. Primary 
pointing, focus and flux calibrations were carried out on Sat- 
urn and secondary pointing calibrations on nearby source 
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Figure 1. Left: Over-density of LBGs in field J1054.7-1254 from tlie ERGS survey (See Douglas et al. 2009). Right: Coverage of 
MAMBO-2 observations, contours display level of coverage in comparison to the best coverage in the field. Shaded area indicates region 
covered by MAMBO-2 observations to ~75% of the best MAMBO-2 coverage (for which 7mm and 12mm observations are in hand). 
LBG positions are marked as circles, with their redshifts, and the dashed box indicates region displayed in Figure 2 



J1058-I-016. Sky dip observations were taken to correct for 
atmospheric opacity and Precipitable Water Vapour (PWV) 
at the site was typically 0.35-0. 7mm (Tau CSO at 225GHz 
~0. 1-0.4). The map was centred on the most densely clus- 
tered region of the field for which 7mm and 12mm observa- 
tions are already in hand (see Davies et al. 2010). A total 
area of ~ 5 arcmin^ was surveyed to 75% of the best total 
coverage encompassing eight LBGs (Figure 1). Data were 
reduced using the dedicated MAMBO-2 reduction package 
MOPSIC (Version 19.01.2010, written by R. Zylka at TRAM) 
utilising an algorithm designed to identify weak sources. In- 
dividual channels are corrected for opacity and instrumen- 
tal noise, and fiux calibrated using data obtained during the 
observations. A Correlated Signal Filter (CSF) was used to 
reduce sky noise, a process designed specifically for the de- 
tection of compact (< 0.5 x the beam size), weak (< 5 x 
rms) sources and which removes diffuse background sources 
of flux. In the resultant map, a root mean square (rms) noise 
of 0.82mJy/beam is obtained over the central ~ 5 arcmin^ 
region. 



3 RESULTS 

Eight spectroscopically-confirmed z ~ 5 LBGs fall within 
the central ~ 5 arcmin^ of the MAMBO-2 observations (Fig- 
ure 2); seven he within a narrow z — 4.95 ±0.08 range while 
the eighth is at z = 5.2. Other spectroscopically-confirmed 
LBGs that fall outside of this central region were excluded 
due to the increased noise in the MAMBO-2 image at larger 
radii (more than twice the rms achieved at the centre of the 
field). Figure 1. No individual LBG is detected at > 2 x rms 
(see figure 2). There are only 4 potential detections at twice 
the rms in this central region, and none at 3x rms, broadly 
consistent with expected statistical fiuctuations. The imme- 
diate fields of each of the four potential detections and the 
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Figure 2. MAMBO-2 map of the field (boxed region figure 1) 
encompassing eight LBGs at 2 ~ 5. Positive flux is indicated as 
light regions. No individual source is detected at a 1.6mjy/beam 
(2 X rms) limit. 



eight LBGs were examined using our optical, near-IR and 
SPITZER/IRAC imaging of the larger fields (Douglas et al. 
2009, 2010) in order to search for optical counterparts and 
assess the potential for source confusion. Only one of the 
four potential detections showed a possible counterpart that 
could account for the signal in the MAMBO-2 data. One of 
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the LBGs had a source within 5" of its position detected 
in all bands from V through to 8 /im which could have con- 
fused the interpretation had that LBG been detected in the 
MAMBO-2 data. 

The typical half-light radii of a UV-luminous stellar 
component of an LBGs in our sample is ~ 0.2" (Douglas 
et al. 2009), while the MAMBO-2 beam has a FWHM of 
10.7" and so covers not just the UV-luminous regions of any 
LBG, but the immediate projected ~ 30kpc around it. It 
is therefore expected that any emission from sources will be 
unresolved in the MAMBO-2 maps even if it is extended 
beyond the UV-emission region, and a 2 x rms upper limit 
on the flux from an individual source can be constrained at 

1.2mm 

< 1.6 mjy 

In order to better constrain the flux from a typical LBG 
and its immediate surroundings (within the beam FWHM), 
data were combined from the positions of the eight LBGs. 
A 35 X 35 pixel (~ 12 x 12 beam size) region centred on 
each LBG is extracted from the map. These 35 x 35 pixel 
images are then combined into an average image taking the 
mean of the eight images at every pixel position. The resul- 
tant image was then convolved with a Gaussian matching 
the MAMBO-2 beam size. The resultant composite image 
is shown in Figure 3. Although there is a region of positive 
flux at the centre of this composite image, it is extended by 
at least two beam widths in one direction, with the peak not 
centred on the position of the objects (beyond the expected 
pointing accuracy of the telescope). Given the limited num- 
ber of independent points in this image (of order 100), the 
significance of this region is not high. We confirm this by 
repeating the stacking analysis using 8 random positions on 
the full MAMBO-2 image. The central ~ 5 arcmin^ region of 
this map contains only ~ 255 independent MAMBO-2 beam 
size regions. If 1000 realisations of a similar combined im- 
age made from random combinations of the 248 positions 
that do not contain a LBG are produced, it is found that 
~ 18% of the realisations contain a > 0.6mJy/beam peak 
in the central position (the flux in the central pixel of Fig- 
ure 3 is 0.6mJy/beam). Consequently, we treat this flux 
as an upper limit to the mean emission of these LBGS, so 
5'i.2mmMEAN < 0.6 mJy, ~ 2.5 times lower than the 2 x RMS 
limit on individual sources. 

The MAMBO-2 beam size (10.7") means that this limit 
not only applies to the FIR flux arising from the rest-frame 
UV-bright regions, but also includes any other emission 
within ~ 5" or ~ 30kpc (and certainly within 1") of the 
sources. This flux may arise from either neighbours at the 
same redshift without signiflcant UV emission, or any under- 
lying system within which individual LBGs are embedded 
as UV-bright star forming regions (see discussion). 

Although the MAMBO-2 field also encompasses the po- 
sitions of several spectroscopically-unconfirmed 2 ~ 5 LBGs 
candidates, we do not include these in our composite anal- 
ysis as if only a small number of these sources are actually 
at 2 ~ 1 — 2 (which can show similar photometric colours 
to z ~ 5 LBGs, see Stanway et al. 2008), they could bias 
an estimate for any FIR flux from true z ~ 5 sources. For 
completeness we also produced a composite image including 
the LBGs which were excluded as having more than twice 
the rms achieved at the centre of the fleld. This compos- 
ite image did not improve upon the depth obtained using 
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Figure 3. Mean composite image of eight LBGs in the field with 
rms'~0.29 mJy/beam. Galaxies are stacked on the optically de- 
rived positions and image is convolved with the MAMBO-2 beam 
size (circle). Contours or 2, 1 and 0.5 sigma are displayed as solid, 
dashed and dotted lines respectively. 



the eight central LBGs - hence we do not discuss it in any 
further analysis. 

3.1 Joint constraints on Lfir, Tdust and Mdust 

In the absence of a detection, we can only jointly constrain 
the FIR luminosity and dust temperature of a typical 2 ~ 5 
LBG, and therefore can only explore a limit on one if we 
make a reasonable assumption about the other. In order to 
relate the limiting flux to a joint constraint, we assume that 
any FIR emission has a grey-body spectrum arising from 
emission by dust heated by starlight. Accordingly we follow 
the procedure used in Aravena et al. (2008) and Stanway 
et al. (2010). We calculate the integrated FIR luminosity 
(Lfir) between 8-1000/im by integrating the grey body spec- 
trum assuming an absorption coefficient of k„ cc ly'^ , as this 
provides the best fit to dust absorption in the spectra of high 
redshift quasars (Priddey & McMahon 2001). We normalise 
this assuming k = 0.4 cm^ g~^ at 250 GHz (See Kruegel & 
Siebenmorgen 1994). In Figure 4, the integrated FIR lumi- 
nosity as a function of dust temperature is displayed for the 
flux limit obtained in the composite image. 

Given the lack of previous study and detections of high 
redshift LBGs in the mm/sub-mm, little is known about the 
appropriate dust temperatures for such sources. Baker et al. 
(2001) use SED fltting of a single, highly lensed LBG at 
z = 2.7 to obtain a dust temperature of T = 33K, and the 
mean dust temperature of submillimeter selected galaxies at 
2 ~ 2 is ~ 35K (Chapman et al. 2005). We note that lower 
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source approaches that of the CMB, the mass required to 
distinguish it from CMB fluctuation increases dramatically. 
In this we assume a source size of 5 kpc, though the re- 
sult does not change significantly for sizes between 2-5kpc 
(typical sizes of 2 ~ 5 LBGs, e.g. Douglas et al. 2009). 

In Figure 5, the dust mass as a function of dust temper- 
ature is displayed for the flux limit obtained in the composite 
image. To constrain the dust mass of a typical source in our 
field, an assumption about the dust temperature must once 
again be made. A dust temperature of 30 K gives a formal 2 
X rms limit to the dust mass of Mdust J; 10* M©. At a higher 
temperatures of 40K, approaching that found for lower red- 
shift ULIRGS, the mass limit is Mdust < 4 x lO'^ Mq. Drop- 
ping the temperature to 20 K leads to a far weaker mass 
limit, Mdust 4 X 10* M©. However this temperature is ap- 
proaching that of the CMB at z ~ 5. Given that there is 
clearly vigorous star-formation in the LBGs, it seems un- 
likely the dust would remain unheated by it or be close to 
equilibrium with the CMB. Hence, the naive expectation is 
that the dust temperature will be higher than this. 



Figure 4. Constants on the Lpijj^ of a typical LBG at 2: ~ 5 
from our stacking analysis of the 1.2mm flux from eight galaxies 
(black); compared to limits derived from 870 /xm observations of 
similar sources (red, Stanway et al. 2010). Solid lines display the 
LpiR a^s a function of temperature for the rms limits. Dashed and 
dotted lines show the 2 X rms and 3 X rms limits. The temperature 
below which dust heating by the CMB occurs is plotted as a dot- 
dashed line. 



redshift (z < 4) ULIRGS are found to have higher temper- 
atures (~ 45 K, Rieke et al. 2009). Although such sources 
have very different redshifts and luminosities to the LBGs 
considered here, they demonstrate that higher temperatures 
are possible. 

In the absence of an exact dust temperature for these 
sources, we assume 30 K as this is close to the temperature 
derived for the 2 — 2.7 LBG - likely to be the most appro- 
priate of the temperatures noted above. Using this, the FIR 
flux limit derived from the composite analysis of a typical 
z — 5 LBG in this work gives a formal 2 x rms limit to the 
FIR luminosity of LpiR < 3 x 10^^ Lq. This limit is several 
times fainter than the FIR luminosities of 2 ~ 3 galaxies 
detected in the sub-mm (Chapman et al. 2005; Chapman & 
Casey 2009; Negrello et al. 2010; Conley et al. 2011), which 
generally have Lfir > 10^^ Lq. With the caveat that the 
temperature is assumed rather than determined, the flux 
limit implies that typical z ~ 5 LBGs have FIR luminosi- 
ties below those of 2 < 4 ULIRGS and intermediate redshift 
BzKs (e.g. see Daddi et al. 2010). 

While assuming a higher temperatures will lead to a 
higher FIR luminosity, it better constrains the associated 
dust mass, which varies inversely with temperature for a 
given flux limit. In order to place a temperature-dependent 
limit on the dust mass, we use the relationship between 
redshift-dependent flux density and dust mass outlined in 
equations 1 and 2 of Weil3 et al. (2007). This includes the ef- 
fect of the increased Cosmic Microwave Background (CMB) 
temperature at high redshift, whereas the temperature of a 



3.2 Obscured star-formation rates 

Star-formation rates (SFRs) in typical 2 ~ 5 LBGs have 
routinely been determined through SED fltting to their rest- 
frame UV-optical emission (sampling star-formation which 
is subject to little dust obscuration), giving typical values 
of few X lOMoyr"^ Verma et al. (2007). While our non- 
detection in the composite image is consistent with a signif- 
icant proportion of star-formation within 2 ~ 5 LBGs being 
relatively free of dust extinction, it is nonetheless interest- 
ing to consider the limit to the dust obscured SFR imposed 
by our non-detection. Using the Kennicutt relation for local 
starburst galaxies (Kennicutt et al. 1998) and our Lfir limit 
we obtain an obscured SFR limit of < 52 Mq yr~^. This sug- 
gests that at least 40% of the total star-formation in 2 ~ 5 
LBGs is unobscured. 



3.3 Combining with LABOCA data 

In addition to investigating the limits derived from this 
study, we can also combine our 1.2 mm observations with 
LABOCA(870 /im) imaging of a similar, 2 ~ 5.16 LBG over- 
density from Stanway et al. (2010) in an attempt to ob- 
tain a limit to or detect the combined emission from the 17 
spectroscopically-conflrmed 2 ~ 5 LBGs across the two data 
sets. We do this by combining the average stacked 1.2mm 
image with a similar average stacked 870 fim image obtained 
by Stanway et al. (2010). We take into account the different 
beam sizes and pixel scales, scaling the 870 /im limit to that 
at 1.2mm by assuming a temperature and spectrum of 30K 
grey body for the average source and accounting for differ- 
ent mean redshifts of the two LBG samples. No source is 
identified at the central position in this combined image at 
a 2 X rms limit of ~0.44 mjy/beam at 1.2mm, representing 
a limit on the average flux of the 17 LBGs. 
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4 DISCUSSION 

4.1 Non-detection in the composite image and 
implications for a ~ 5 sources 

4.1.1 Model predictions for FIR emission ~ 5 LBGs 

Until recently, there have been no self-consistent model- 
based predictions for the expected fluxes of sources such 
as these. After these observations were carried out, Gon- 
zalez et al. (2011) made predictions for the 850/im flux of 
distant LBGs as a function of their UV luminosities, based 
on the GALFORM semi-analytic model of galaxy evolution. 
Although they concentrate on redshifts « = 3 and z = Q, 
the small variation in prediction at around at rest-frame 
UV absolute magnitude, Mi5()o(^B)^51og(h) ~ —19.8 (the 
typical value for the ~ 1500A absolute magnitude of the 
ERGS sources, where h is the dimcnsionless Hubble param- 
eter) indicates we can use these predictions at z ~ 5. Look- 
ing at their figure 11, the median flux predicted is around 
0.4 mjy at 850 /im, which we calculate to bo equivalent to 
0.2- 0.25 mjy at 1.2 mm (for grey body temperatures in the 
range 30-50K). 

We also infer from their figure 11 that about ten per 
cent of their sources should be brighter than ~ 1 mJy at 
850 ^m. The limits determined in both this and the earlier 
work of Stanway et al. (2010) for the typical source, are 
roughly double the median value predicted by Gonzalez et 
al. (2011). While our non-detections are compatible with 
the model, it is clear from our attempt at combining both 
the current and Stanway et al. (2010) data sets that a sam- 
ple size double that of the current one (observed to similar 
depths or deeper) will start to challenge or validate these 
model predictions. Of course, this is just one model with 
one set of assumptions and prescriptions; having the abil- 
ity to test more models which use different methodologies 
would increase the utility of such observations. Interestingly 
an alternate model by Finkelstein et al. (2009), based on the 
assumption that sub-mm emission in these sources is pow- 
ered by the UV-luminous component, gives predictions that 
span a range covering both the predictions of Gonzalez et 
al. (2011) and our limits on the typical LBG. 

4.1.2 The nature and baryonic content 0/ 2: ~ 5 LBGs 

As noted earlier, we obviously cannot know the dust temper- 
ature for any material within these galaxies in the absence 
of detections at multiple frequencies. However, by assuming 
a temperature we are able to explore how the limit to the 
dust mass relates to the mass in stars for these sources. For 
temperatures of 30 K or above, a dust mass limit for LBGs 
at ^: ~ 5 is found to be roughly 1% to 10% of the typi- 
cal mass of their stellar component (depending on study, 
for a typical 2 ~ 5 LBG similar to those observed here, 
M* ~ few X 10^"^° Mq at least, Verma et al. 2007; McLure 
et al. 2009). In addition, the molecular gas content of a 
typical LBGs in this field has been found to be compa- 
rable to this typical stellar mass, at most (Mh2 < a few 
xlO^ M0, assuming a CO flux to H2 mass conversion 'X- 
factor' = 0.8 M0(K km s~^ pc^)^"*^ for consistency with ear- 
lier work, see Solomon & Vanden Bout 2005; Davies et al. 
2010, and references therein). Hence, assuming a dust tem- 
perature of > 30K, no more than ~ 50% of the total baryonic 



mass of LBGs at z ~ 5 is in their dust and molecular gas 
component - to the current limits and given the assumptions 
above. Clearly, using a larger CO flux to H2 mass conver- 
sion factor will increase the upper limit to the molecular gas 
content. Using an X-factor more typical of low redshift spi- 
rals, 4.5M0 (Kkms~^ pc^)~^, with these limits, allows up 
to ~ 80% of the baryonic mass of 2; ~ 5 LBGs to be in the 
form of molecular gas and dust if a typical source has a stel- 
lar mass of ~ IO'^Mq. For sources with stellar masses closer 
to IO^^Mq, even with this X-factor, the gas would account 
for less than 50% of the baryonic mass. As we know little of 
the correct CO flux to H2 mass conversion factor for these 
sources we can not constrain this fraction further and must 
consider this limitation in any conclusions drawn from such 
an analysis. 

By determining limits to the total baryonic content of 
2 ~ 5 LBGs we can explore the nature of the LBG phe- 
nomenon at those redshifts. A substantial fraction of these 
objects show distorted or multiple UV components on sub- 
arcsecond scales, corresponding to a few kpc (e.g. Conselice 
& Arnold 2009; Douglas et al. 2010). There has been detailed 
discussion in the literature as to what these morphologies 
tell us about the nature of LBGs. Through a comparison 
to low redshift LBG analogues, Overzier et al. (2009) show 
that such morphologies can be explained as the result of gas 
rich mergers of low mass systems triggering starbursts, while 
Law et al. (2007) explore the possibility that the UV compo- 
nents are relatively unobscured regions of star-formation in 
a larger underlying gala:xy which is either intrinsically less 
UV-luminous or is extincted. 

Given the separation between components (typically 
~ 0.5", a few kpc), any encompassing or underlying galaxies 
clearly have to bo larger and more massive than the sizes and 
masses attributed to the UV-luminous regions (so with a to- 
tal stellar mass > IO^^Mq and possibly considerably larger). 
Both a gas rich merger or a larger xmderlying system could 
easily lead to appreciable far infra-red (FIR) luminosities, 
through significant rapid dust production in the case of the 
merger or simply by the amount of obscuring material/total 
mass in the case of the underlying system. However, it is 
perhaps easier to explain a lower FIR luminosity in the case 
of a merger, as it could be argued that in these LBGs the UV 
emission traces the majority of the material in the system. 

The apparent lack of large quantities of material that 
could be detected in the millimetre and sub-millimetre (dust 
and molecular gas) in these SA'stcms weakens the case for 
LBGs being unobscured super starburst regions embedded 
in much larger systems, (see discussion in Davies et al. 2010). 
While this also argues against very gas-rich mergers that 
promptly produce large amounts of dust, it poses no con- 
straint on mergers where most of the baryonic material is 
traced by the UV emission and is in the form of stars. 

4.2 Comparisons to lower redshift populations 
and implications for future observations 

Given the lack of detections in the current data it is worth- 
while exploring how much deeper further observations may 
need to probe before the population is detected in emission 
from dust. We can do this by investigating the properties of 
lower redshift {z < 4) sources which, while not ideal, may 
have use as analogues for 2 ~ 5 LBGs. 
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Figure 5. Constants on the dust mass of a typical LBG at z ~ 5 
from our stacking analysis of tlie 1.2mm flux from eight galaxies 
(black) ; compared to limits derived from 870 /xm observations of 
similar sources (red, Stanway et al. 2010). Solid lines display 
the dust mass as a function of temperature for the rms limits. 
Dashed and dotted lines show the 2 X rms and 3 X rms limits. 
The temperature below which dust heating by the CMB occurs 
is plotted as a dot-dashed line. 



The limit on the stellar mass to cold dust ratio ob- 
tained here (> 10) is similar those observed in some local 
starbursts. In a study of metal poor blue compact dwarf 
galaxies, Hunt et al. (2005) derive stellar to dust mass ra- 
tios of 7.5 and 13 for two sources. Both stellar and dust 
masses in Hunt et al. (2005) are derived from SED fitting 
and no errors are given, hence these are only best-fit val- 
ues. Although these sources are less massive than LBGs at 
z ~ 5, they have comparable specific star- format ion rates 
and metallicity (0.14 Zq and 0.2 Z© compared to 0.1 - 0.2 Z0 
obtained for a ~ 5 LBGs from SED fitting, e.g Verma et al. 
2007, and rest-frame UV spectral slope analysis, Douglas et 
al. 2010). In addition, these local starbursts have emission 
dominated by small-grain, Type-II SNe formed dust as is 
present in high-z sources (Maiolino et al. 2004) and there- 
fore, they should display similar dust characteristics. 

If these systems are good models for 2: ~ 5 LBGs, obser- 
vations not much deeper than those presented here should 
start to detect the typical LBG at FIR wavelengths. For 
example, using a stellar mass of M = lO" Mq and dust 
temperature of > 30K, a stellar to dust mass ratio of < 13 
predicts a dust mass of Mdust 8 x 10^ M©. Assuming the 
same grey body model outlined above, an rms of ~ 0.14 mjy 
is required for a 2 cr detection of such a source. For a similar 
MAMBO-2 observation and stacking the 8 LBG in the field 
a ~ 30h of integration time are required to reach this limit. 
In comparison, for a similar integration time to that in this 
study, a composite image of > 30 LBGs would be required 
to obtain a detection. 

A more direct comparison can be made with UV- 
selected galaxies at lower redshift. LBGs at z ~ 3 are iden- 



tified via bright rest-frame UV continuum emission using a 
similar method to those at z ~ 5. Their UV/optical proper- 
ties indicate they have comparable star-formation rates, but 
with higher metallicities, slightly larger dust extinctions in 
the rest-frame UV and have typically older and more mas- 
sive detectable stellar population (by close to an order of 
magnitude) than z ~ 5 LBGs (see Verma et al. 2007), so 
their use as analogues is limited. 

Rigopoulou et al. (2010) detect the average emission 
from IR-luminous z ~ 3 LBGs at 250/^m with Herschel- 
SPIRE and calculate Lfir = 2.8 X 10^2 L© for Td^st = 45K. 
Looking at figure 4, we note that if we assumed this tem- 
perature, our flux limit already probes this luminosity at 
z ~ 5. However, for sources not selected to be IR-luminous, 
the picture is different. Thermal dust emission has been de- 
tected in several strongly lensed LBGs at z ~ 3 giving typ- 
ical dust masses of ~ few x lO^M© (e.g. Baker et al. 2001; 
Coppin et al. 2007). Stanway et al. (2010) calculate a sim- 
ilar limit to the dust mass of typical z ~ 3 LBGs from the 
flux limit obtained by Webb et al. (2003) from observations 
of such sources in the Canada-UK Deep Submillimeter Sur- 
vey (again assuming Tdust = 30K). If these values are also 
appropriate at higher redshift, then significantly deeper ob- 
servations will be required to detect the dust emission from 
LBGs at z ~ 5, observations which may run into the confu- 
sion limit of single-dish observations. 

However, the required depth will be straightforwardly 
achievable with ALMA when it reaches full science mode. 
Figure 6 displays the dust mass limit as a function of inte- 
gration time for ALMA in full science mode assuming a grey 
body at 30 K. These values are calculated using the ALMA 
exposure time calculator observing a field at a declination of 
-12 deg at 250 GHz, with 7.5 GHz bandwidth, T,fcj,=44.3K 
and Tsys = 132.5K. Clearly ALMA will make the detection 
of FIR emission from such sources routine, obtaining the 
same limit as our MAMBO-2 observations in < 10 seconds 
and detecting Mdust ~ 10^ M© sources at z~ 5 in just tens of 
minutes. While this represents a relatively short exposure for 
ALMA, it does mean that these sources will have to be tar- 
geted individually rather than as part of anything other than 
a deep blind survey. The source density of spectroscopically- 
confirmed ERGS LBGs is approximately one per 6 arcmin'^, 
or one source per ~ 50 ALMA primary beams at 250 GHz. 
While fainter LBGs are more numerous, they are also likely 
to be fainter in the mm/sub-mm {e.g. Gonzalez et al. 2011). 



4.3 Overlap between LBG and SMG populations 

at z ~ 5 

Several spectroscopically-confirmed z > 4 SMGs have been 
identified in the past three years {e.g. see the list in Kohno 
2011). While the optical photometry of these sources is in- 
complete in the literature, of the seven sources in Kohno 
(2011), five would have been selected in a survey using the 
Lyman break technique (see Schinnerer et al. 2008; Coppin 
et al. 2009; Daddi et al. 2009a,b; Riechers et al. 2010; Knud- 
sen et al. 2010). Given that a large fraction of distant SMGs 
are potentially detectable in LBG surveys (even if for some 
of them the UV-light is contaminated by AGN emission, un- 
like most LBGs, e.g. Coppin et al. 2009), it is informative to 
quantify the overlap between the population of typical z ~ 5 
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Figure 6. The detection limits obtainable as a function of inte- 
gration time for the full sensitivity (50 antennae) ALMA array. 
The ALMA exposure time calculator was used to determine the 
2xrms for a range of exposure times. This was converted to a 
dust mass limit assuming a 30 K grey body as in the rest of this 
work. The solid diagonal line displays limits achievable for a sin- 
gle source and dotted line the limits achievable when stacking 8 
sources, as carried out here. The red vertical dashed line indicates 
the mass limit derived from these MAMBO-2 observations. The 
blue dot-dashed line displays the predicted dust mass if 2 ~ 5 
LBGs have the same M* /M^dust blue compact dwarf galaxies 
(Hunt et al. 2005). The orange dashed-dot-dot line displays the 
Gonzalez et al. (2011) predictions for 2 ~ 5 LBGs. The green 
large dashed line displays the dust mass limit calculated from the 
Webb et al. (2003) sample of z ~ 3 LBGs. The black vertical 
dashed line is a predicted lower limit determined by applying the 
Kennicutt relation to a star-formation rate of ~ lOM0/yr, typi- 
cal for the UV-detected star-formation rates in z ~ 5 LBGs, and 
converting to a dust mass assuming a Tdust=30¥i. 



LBGs and the more extreme (as far as star-formation rates 
are concerned) SMGs at this redshift. 

We have shown that the typical 2 ~ 5 LBG is not a 
strong FIR source. By combining our results with those of 
Stanway et al. (2010) we can place a limit on the fraction 
of LBGs that could be mJy sources at ~ 1 mm. None of 
the 17 spectroscopically-confirmed sources probed are de- 
tected individually at 1-3 mJy and together have an aver- 
age flux of less than 1 mJy implying that less than 10% of 
these confirmed LBGs have submillimeter flux at the milli- 
Jansky level. In addition to these sources, both clustered 
fields contain multiple photometrically-selected LBG can- 
didates without spectroscopic confirmation, either because 
they were not targeted by spectroscopy due to observational 
constraints, or because spectroscopy was inconclusive (Dou- 
glas et al. 2009, 2010). None of these were detected at 870 /im 
or 1.2 mm. Accounting for the fraction of these candidates 
that will be at z ~ 5 (Douglas et al. 2009), this almost dou- 
bles the number of 2 ~ 5 LBGs covered in the deep regions 
of the two fields, indicating that less than 5 % of 2; ~ 5 LBGs 
have flux at 2 - 3 mJy. 



In order to place an upper limit on the fraction of FIR 
bright sources in the 2 ~ 5 LBG population, we can ex- 
trapolate further. Given the volume probed and the lack of 
detections in both this and the Stanway et al. (2010) work, 
if the typical source is detected at just below our limits (this 
does not have to be true but will maximise the fraction of 
bright sources) and if at this flux level we probe the expo- 
nential (bright) end of the LBG FIR luminosity function, 
we expect less than 1 % of the spectroscopically confirmed 
2 ~ 5 LBG population to give a > 5 mJy detection. If the 
typical source is considerably fainter, then an even smaller 
fraction will be this bright. 

Furthermore, it should be remembered that these 
sources are selected from the two most over-dense fields of 
the ten similarly sized ERGS pointings. If there is any bias 
towards more massive starburst galaxies being preferentially 
found in the regions with the most evolved large scale struc- 
ture at high redshift (the most massive systems at any epoch 
are likely to be found in the most over-dense environments, 
e.g. see Mo & White 2002; Blain et al. 2004), the lack of a 
detection in these fields may indicate an even lower limit on 
the fraction of potential FIR bright LBGs in this and in any 
unbiased LBG survey as a whole. 

4.4 The lack of other detections in the field 

Given the expected typical sub-mm flux of LBGs and known 
sub-mm number counts, it is perhaps unsurprising that no 
sub-mm bright sources were detected within this fleld to 
our limiting flux (at ~ 700 sources/deg^, we would expect 
~ 1 ± 1 source in a randomly chosen field this size, based 
on the numbers in Hatsukade et al. 2011). However, this 
field was selected because of its unusual clustering of 2 ~ 5 
LBGs, indicating that it contains a significantly over-dense, 
and therefore evolved, structure at this redshift. 

It has been argued that luminous SMGs require signif- 
icant early hierarchical evolution of their environments to 
exist and hence should be found in high density peaks of 
the matter distribution at high redshift (e.g. Blain et al. 
2004). Recently, a 2 ~ 5.3 SMG was discovered within an 
over-density of photometrically-selected LBGs at the same 
redshift (Capak et al. 2011), supporting this picture. The 
LBG over-densities in the currently studied field and in that 
studied by Stanway et al. (2010) are comparably dense. All 
three fields could potentially be the progenitors of low red- 
shift massive clusters, or at least go on to form very massive 
galaxies by the present-day. Therefore, they are equally plau- 
sible regions in which to search for massive (UV-obscured) 
submillimeter galaxies. 

Despite the LBG clustering indicating a high density 
environment, we do not detect any sources within either 
field with fiuxes above a few mJy (with FIR luminosities 
of more than a few xIO^^Lq and therefore, dust masses 
of above a few xlO^M© assuming T = 30K, and inferred 
star formation rates above a few x 100 M© yr~^). A possible 
explanation for this is that the FIR luminous stage giving 
rise to SMGs is expected to be short (~100Myr). Therefore, 
it is entirely possible that at least one such system, which is 
at an evolutionary stage when it is not FIR luminous, will 
be in the extended structure rendering it undetectable in 
our observations. The typical timescale of the LBG phase at 
these redshifts is comparable or less (< 100 Myr, Verma et al. 
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2007; McLure et al. 2009). However, as discussed in Stanway 
et al. (2008) and Douglas et al. (2010), given the relatively 
high volume density of such sources, as some LBGs fade 
in the UV, others "switch on". This allows the large-scale 
structure to bo continually traced by LBGs, albeit different 
LBGs at different times. Given the lower volume density 
of SMGs it is possible that we have simply observed the 
structure at a time when none of the possible SMG systems 
are FIR luminous. 

These MAMBO-2 observations have not probed the en- 
tire spatial extent of the over-density identified in the optical 
data (see Figure 1) and it is possible that an SMG similar to 
that identified in Capak ct al. (2011) could lie beyond the 
high sensitivity region in Figure 1. However, the observa- 
tions are centred on the most over-dense region in the field, 
perhaps the most likely position for such an object. Obvi- 
ously the sensitivity of these observations is such that an 
object with a FIR luminosity (and star-formation rate, dust 
mass etc) a few times lower than the known z ~ 5 SMGs 
would not have been detected. Given that the known objects 
represent the extreme of the population, we cannot rule out 
this possibility. 



5 CONCLUSIONS 

Constraining the thermal dust component of LBGs at z ~ 
5 is an important step in developing a picture of star- 
formation and galaxy evolution at high redshift. We have 
carried out 1.2mm MAMBO-2 observations of a field over- 
dense in LBGs at « ~ 5, allowing the simultaneous obser- 
vation of eight LBGs. No individual source is detected at 
a 1.6mJy/beam (2 x rms) limit. When stacking the flux 
from the positions of all eight galaxies we obtain a limit of 
0.6 mjy/beam for the average emission from the objects. As- 
suming a dust temperature of T = 30 K this corresponds to 
a FIR luminosity and dust mass limit for a typical LBG at 
« ~ 5 of Lfir < 3 X 10" Lq and Mdust < 10* M©, less than 
ten per cent of the stellar mass (~ 10^ M©). However, this 
dust mass limit is dependant on dust temperature, which is 
yet to be constrained in high redshift galaxies. We estimate a 
limit to the total star-formation rate in a typical « ~ 5 LBG 
of < 52M0 yr~^ using the Kennicutt relation. This compares 
to minimum star-formation rates of a few xlOM0yr~^ de- 
termined from the rest-frame UV emission of typical z ^ 5 
LBGs. When combined with our previous results on a simi- 
larly clustered field (Stanway et al. 2010), we have observed 
17 spectroscopically-confirmed ^ ~ 5 LBGs with uniform re- 
sults, ensuring that these studies robustly characterise this 
population. 

No other objects are identified in the field, ruling out 
any source similar to the known « ~ 5 SMGs lying in the 
same, highly over-dense structure as the observed LBGs. 

The limits on the dust mass within the LBGs and their 
immediate (~ 30kpc) environments weakens the possibility 
that these objects arc relatively unobscured super starburst 
regions embedded in dust-obscured more massive, larger sys- 
tems, especially when combined with the limit on the molec- 
ular gas emission from our earlier work. 

Between this and our previous work we have charac- 
terised the molecular gas and dust properties of the typical 
population of 2: ~ 5 LBGs to the limit of the available instru- 



mentation. With the advent of ALMA we should be able to 
straightforwardly detect both the cool gas and dust phases 
of these LBGs and their environments (in targeted observa- 
tions) , even if they prove to have a baryon content towards 
the lower end of the predicted expectations. 
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